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We present the calculation of the fermionic contribution to the QCD two-loop amplitude for e + e~ — > qqg. 



1. Introduction 

Today's high energy collider experiments have 
reached such an accuracy that next-to-leading order 
(NLO) theoretical predictions are no longer sufficient. 
While for inclusive observables like for example the 
total hadronic cross section in e + e~ annihilation the 
step from NLO to next-to-next-leading order (NNLO) 
is far less demanding - and has been taken a long time 
ago - the contrary is true for less inclusive quanti- 
ties as for example jet rates. Given these difficulties 
it is clear that only for specific reactions NNLO cal- 
culations can be envisaged. Important reactions for 
which one should go beyond the NLO approximation 
are for example: Bhabha scattering, pp,pp — > 2 jets, 
and e + e~ — > 3 jets. Among these examples the e + e~ 
annihilation into 3 jets is of particular interest. His- 
torically this reaction was one of the first clean tests 
of Quantum Chromodynamics (QCD) as the theory 
of strong interaction. Today the 3-jet production in 
e + e~ annihilation is still an excellent laboratory for 
precise tests of QCD and jet physics. If one stud- 
ies how a jet is modeled in the perturbative predic- 
tion, it also becomes obvious that only in higher or- 
der perturbation theory will give a reliable descrip- 
tion (c.f. fig. [j]). Due to the fact that lepton collid- 
ers provide a very clean environment for QCD anal- 
yses e + e~ — > 3 jets is also an excellent reaction for 
measureing the QCD coupling oc v with high accuracy. 
It is exactly here where the need for NNLO predic- 
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Figure 1 . Modelling of a jet in leading a), next-to- 
leading b), and next-to-next-to-leading order c). 



tions for 3-jet production becomes most prominent: 
at present the measurements of oc s are plagued by 
scale uncertainties of the theoretical predictions ^ 
(see also H. Stenzel, these proceedings). These un- 
certainties arise from uncalculated higher order con- 
tributions. They can be reduced by a NNLO calcula- 
tion. For a 3-jet prediction at NNLO accuracy several 
ingredients are needed: First the matrix elements for 
e + e~ — > qqggg, qqq'q'g nave to be evaluated. They 
are obtained from a leading-order calculation and are 
known for a long time [g, [|]. Second the one-loop 
matrix elements for e + e~ — > qqgg,qqq' q' are needed. 
They have been calculated in refs. [[j], || [1 Hi ■ Finally 
one needs the NNLO amplitudes for e + e~ — > qqg 
[[§|,0[ll]]. In the following the calculation of the 
NNLO matrix elements is reported. One should keep 
in mind that the combination of these 3 ingredients 
into an infrared finite cross section is still a highly 
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non-trivial - and so far unsolved - problem. 

2. Old and new approaches to calculate NNLO 
amplitudes 

Due to the tremendous activities in the field enor- 
mous progress has been made in the past. All the 
relevant double box scalar integrals have been calcu- 
lated [12 13|,[l4],[l5[]. For a few reactions these results 
have been used to obtain NNLO scattering amplitudes 

Today one 

can say that a standard approach to perform these cal- 
culations exist. Essentially it consists of three steps. 
First one generates automatically the relevant Feyn- 
man diagrams. In a second step the reduction of ten- 
sor integrals to scalar integrals is done via Schwinger 
parametrization [ gl|g7[ ], The scalar integrals obtained 
in this step appear in shifted dimensions and with 
raised powers of the propagators. In a last step, equa- 
tions derived from partial integration and Lorentz in- 
variance are used to reduce the scalar integrals to a 
small set of master integrals [ |28| , |30| ] . The final 
result is than obtained in terms of these master inte- 
grals. They have to be calculated analytically. The 
bottleneck of this approach is that one has to create 
and to solve a huge system of equations to obtain the 
desired reduction to master integrals. In contrast to 
the naive expectation, topologies which at first glance 
look not so difficult may involve much more work 
than the more complicated once like for example the 
double box. In ref. [^TJ we have proposed a different 
method. The basic idea is that one applies only obvi- 
ous reductions like for example the triangle rule. In 
particular one does not perform a complete reduction 
to master integrals. Instead one calculates the scalar 
integrals in higher dimensions and with raised powers 
of the propagators directly. To illustrate the method 
let us consider the so-called penta-box. A correspond- 
ing Feynman diagram is shown in fig. ||a. Applying 
Schwinger parametrization one arrives at the scalar 
topology shown schematically in fig. ||b. An obvious 
application of the triangle rule allows the elimination 
of two propagators [p2[], The resulting topology is 
often called C-topology. This reduction is shown in 
fig. ||. As mentioned earlier the C-topology appears 
in higher dimensions [d = 2m — 2e, m 6 N) and with 
raised powers of the propagators (vi , . . . ,Vs). In the 
approach advocated here an analytic expression for 
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Figure 2. Tensor reduction of the general pentabox 
topology a) to scalar integrals in higher dimensions 
and raised powers of the propagators b). 




Figure 3. Application of the triangle rule yielding the 
C-topology C(w = & + e, Vi , v 2 , v 3 , v 4 , v 5 ) 



the C-topology is needed. Such an expression can be 
obtained in terms of infinite sums [f3l|]. The generic 
structure is shown in eq. ([!]). 



C(m=- + e,v 1 ,...,v 5 )~ £ £ ^-^ 



l[=Ol 2 =0 



2m-2£-v 12 35 2m-2e-V 2 345 + 2m — 2e — V125) 



1 2 r(ii + l+2m-2e-Vi235) 

rp'i +h + 3m- 3e-v 12345) 
X r(r 2 + l + 2m-2e-V2345) 

rp'i + 2 2 + 2m - 2e - V 2 3s)r(i2 + 2m - 2e - V 345 ) 
r(u+«2 + 4m-4e-Vi2345-V5) 

(1) 

It is obvious that such a representation is completely 
useless unless one is able to evaluate the e-expansion 
in terms of known functions. Fortunately this is possi- 
ble using a few basic algorithms for nested sums [pljl. 
The basic idea here is that one rewrites the sums using 
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the relation 



n-l 



52 52 fl y = «oo + 52 a °/ + L a '° + L L a i{n-jyi 2 ) 

i'=0 7=0 ;'=1 i'=l n=l;'=l 

Furthermore one can expand the T-functions by the 
means of 



r(-n + e) 



-1)" 



en! 



r(l+e)(l+Si(n)e 



+Si,i(n)t? + S u ,i(n)e 3 + S u ,i,i(n)e 4 + ... (3) 

and related identities [[33"|]. Here Si ,...(«) denote the 
harmonic sums, see e.g. ref. j33|]. Inspecting the sums 
which appear using eq. (|j) and eq. (Q) one observes 
that the algorithms needed are just a generalization 
of the algorithms studied in ref. p3|]. In particular 
defining the nested sums as a generalization of the 
harmonic sums by 

11 1*1 >k-l -Jl x 'k 

S(n;m h ...,m k ;xi,...,x k ) = 52 52 ' " L W ' " W 

i'l = li 2 = l i fc =l 'l 'it 



1=1 



(4) 



the following four basic operations are sufficient to 
evaluate the e-expansion of the C-topology in terms 
of multiple polylogarithms p4]: 



1. Multiplication 

S(n;m\,...;x\,...) x S(n;m[, ...;x[, ...) (5) 

2. Sums involving i and « — i: 

«-l j.i 

£-±-S(i;/n 2 ...;*2,...) 

i=i ' 

-5(n — i;m 2 , ■ ■■',x / 2 , ...) (6) 
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3. Conjugations: 

17 / \ 

% y i j(- i y^ s ( i ' m i>"-> m k> x u---,Xk) (7) 

4. Sums involving binomials and n — i: 

n-l 



(8) 



£^J(-iy^%m 2 ..,^...) 



(n -0' 



— m' 2 , ...;x 2 , ...) 



In ref. pl| ] we have given explicit algorithms for 
the four basic operations. These algorithms are well 
suited for an implementation in a computer algebra 
system. For their implementation we used Form3 [ |35| | 
and Ginac |36|, ^7||. Using these programs we cal- 
culated more than 300 C-topologies analytically. In 
particular we have checked that our results for the 
two master integrals C(2, 1,1,1, 1,1), C(2, 1,1, 1,1,2) 
agree with those given by Gehrmann and Remiddi 
||38||. All the simpler topologies can be calculated in 
the same way. 

3. Results 



The amplitude Ay for the process 
f e~->-Y* <l(pl)g(p2)q(p3) 



(9) 



can be written as a leptonic part multiplied by a 
hadronic part H^. 



Jl V 
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= i ft (Hf + g^ + (g)V , + ...W(lQ) 

with s = (pi + p2 + pi) 2 and g s = ^/4na s . Further- 
more the hadronic part can be decomposed according 
to the colour structure. In particular we have 

rr(2) _ T a ( N 2rj(Tj , „(2) , 1 „(2) 



N 



(11) 



with T"g the generator of the SU (N) gauge group and 
«/ the number of massless quarks. The last contri- 
bution arise from a closed quark loop coupled di- 
rectly to the photon. (The factor Hf involves a sum 
over the charges of the active quarks.) We note that 
the contribution proportional to rth is completly de- 
termined by the renormalization. It does not involve 
any non-trivial two-loop topology. In this talk we dis- 
cuss only the contributions proportional to n/N and 
-*£. The hadronic contribution can also be decom- 
posed according to the spinor structure: 



Cl-(Pl\h\Pl) E ; 
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+C2-^{pl\fa\P3)(e g ' Pl)P3p 
+ c 3-^(Pl\h\P3)(£g- P3)Pl M 
+C 4 ^(pi\p'2\p?,)(Zg- Pl)pi M 
+ --- + Cl3-^(Pl\h\P3)(£g- P3)P2v (12) 

where the scalar functions ci — C13 depend only on the 
ratios x\ = 2(p\ ■ pz)/s, x 2 = 2{j> 2 ■ Pi)/s. Due to vari- 
ous constraints arising from gauge invariance or quark 
anti-quark antisymmetry only 4 of them are needed. 
All the remaining ones can be expressed in terms of 
these 4 functions which we chose to be C2 , C4 , and 
c\2- The perturbative expansion in a s of the functions 
c, is defined through 



g .(c?> + (±)c?h(%) c?h0{al) .(13) 



Using the methods described in the previous section 
we calculated all the 13 functions and used the con- 
straints as an additional check of our calculation. Fur- 
thermore we checked the structure of the ultraviolet 
as well as the infrared divergencies [39]. We have 



also compared our results with results published re- 
cently Inrefs. gg] the result is written in terms 
of 2-dimensional harmonic polylogarithms instead of 
the multiple polylogarithms used here. Expressing 
the 2-dimensional harmonic polylogarithms in terms 
of multiple polylogarithms we found complete agree- 
ment for the n fN and ^- contribution calculated here. 

Following ref. [g] we used the general structure of 
the infrared divergencies as described by Catani 
to define the finite parts of the coefficients cf. 



JZ)fin _ (l),ren 



_l(i)( e ) c «. ren _l(2)(e) 



„(0) 



(14) 



with the one- and two-loop insertion operators iW (e) 
and l( 2 )(e) given in ref. [p9|]. As an example, we 
present the result for the n /A^-contribution to the fi- 
nite part 0?^° at two loops, 



„(2),fin 



(xi,x 2 ) 
ln(x 2 ) 



_ / 19 1 
N ~ Hf V 18*i(l-JC2) 



36 (1-Jfc) 



(Xi + X 2 ) 



59 



50 



x\{\ -x 2 ) 



-12 



(l—3xi)(l— jci) 



xi(l -X1—X2) 



(xi-x 2 ) 
(x\ + x 2 ) 3 



Ri(xi,x 2 ) 



5 (2-x 2 ) | 2 
12xi(l — X2) 2 (xi+x 2 ) 2 



Li 2 (l -x 2 ) 



1 {2-x 2 ) 



ln(* 2 ) 2 + -C 2 



4x^1 -x 2 ) 2 
1 (7xi — 5x 2 — 6x\X2 + 5x2 2 +x\ 2 ) 
12 (1 -x 2 )x\(l —Xi -x 2 )(xi +x 2 ) 



ln(xi) 



R(xi,x 2 ) 
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6 XZ - 12 



(1 -xi -x 2 )(l -x 2 )(x\ +x 2 
19 1 (1+xQ 1 (l-3x 1 )(l-x 1 ) 
12 12(1 -x 2 ) 3 (l-xi-x 2 ) 
X[(l +x\) 



+ 5xl_5_x 2+2 

12xi 6 xi (xi +x 2 ) 
1 (2-x 2 ) 2 



12xi(l-x 2 ) 2 (xj+x 2 ) 2 
(2-x 2 )ln(x 2 ) 1 



Li 2 (l -x\] 
1 



(15) 



2xi(l— x 2 ) 2 2xi(l— x 2 ) 

We have introduced the function R(xi,x 2 ), which is 
well known from [5oj], 

R(xj,x 2 ) = Qln(x 1 )ln(x 2 )-ln(x 1 )ln(l-xi) 



+ -C 2 -Li 2 (xi) ) +(x\ 



x 2 ) 



(16) 



In addition, it is convenient, to define the symmetric 
function Ri (xi ,x 2 ), which contains a particular com- 
bination of multiple polylogarithms [B4J, 



Ri(xi,x 2 ) = ln(xi)Li 



U.i 



X\ 



X\ +X2 ' 



-,xi+x 2 



-~£ 2 ln(l -xi -x 2 ) - ~ ln(xi) ln(x 2 ) ln(l -xi -x 2 ) 



-ln(xi)Li 2 (xi+x 2 ) — Lii.2 



x\ 



xi+x 2 



,xi+x 2 



A'l 



Li 2 ,i 

\x\-\-x 2 
+ (xi <^x 2 ) 



,x\+x 2 +Li3(xi+x 2 



(17) 



As one can see all multiple polylogarithms have sim- 
ple arguments. As a consequence it is straightforward 
to obtain the analytic continuation, which can be used 
for the crossing of the amplitude. 

4. Conclusion 

In this talk we have presented the calculation of a 
specific colour structure of the QCD two-loop ampli- 
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tude e + e~ — » ggg. The computation has been done 
by a completely new method developed by the au- 
thors [f31[]. It has been shown that the approach which 
is based on nested sums and their algebraic proper- 
ties provides a very efficient method for two-loop 
calculations. Furthermore the tools developed in this 
approach can also be used for a systematic expansion 
of generalized hypergeometric functions. In addition 
the presented results provide a very important cross 
check on recently obtained results 
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